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ABSTRACT 

We investigate what we can learn about galaxy formation and evolution from the data which will be ob- 
tained by the forthcoming large submillimeter/millimeter facilities, mainly by the Atacama Submillimeter Tele- 
scope Experiment (ASTE) and the Atacama Large Millimeter Array/Large Millimeter and Submillimeter Array 
(ALMA/LMSA). We first calculate the source counts from 350 jim to 3 mm using the empirical infrared galaxy 
number count model of Takeuchi et al. (2001). Based on the number counts, we evaluate the source confusion and 
determine the confusion limit at various wavebands as a function of the characteristic beam size. At submillime- 
ter wavelengths, source confusion with the 10 - 15-m class facilities becomes severe at 0.1 to 1 mJy level, and 
astrometry and flux measurements will be difficult. However, we show that very a large-area survey of submil- 
limeter sources brighter than 10 - 50 mJy can provide a unique constraint on infrared galaxy evolution at z = 1-2, 
and such a survey is suitable for the ASTE. In addition, such a survey enables us to study the clustering properties 
of submillimeter sources, which is still highly unknown. We also find that the 5er-confusion limit of the LMSA 
is fainter than 1 /iJy, which enables us to study the contribution of sources at extremely large redshift. When 
we discuss such a deep flux limit, the dynamic range of a detector should be taken into account, since extremely 
bright sources make it impossible to detect the faintest sources near the detection limit. We evaluate the probabil- 
ity that sources that are 10 3 times brighter than the 5c-detection limit of the LMSA and the ALMA exist in the 
field of view. We find that the probability is < 3 x 10 -4 , and therefore we do not have to worry about the dynamic 
range. The source counts at such faint flux levels give important information of the epoch of galaxy formation. 
We then show that multiband photometry from the infrared (by ASTRO-F) to the millimeter can be utilized as a 
redshift estimator. We examined the performance of this method by Monte Carlo simulations and found that it 
successfully works if we have reasonable measurement accuracy. In addition, we compare the observed 1.4, 5, 
and 8-GHz source counts with our model counts to examine the contribution of star forming galaxies to the faint 
radio galaxies. We find that the faintest radio number counts (~ 1 /iJy) are dominated by actively star-forming 
galaxies which lie at intermediate redshift z ~ 1-2. 

Subject headings: galaxies: evolution — galaxies: formation — infrared: galaxies — radio continuum: galaxies 
— submillimeter 



1. INTRODUCTION 

A significant fraction of galaxies in their youth are expected 
to be dusty starburst since the first epoch of star formation pro- 
duces a lot of dust. The radiation of underlying stars is absorbed 
by dust that re-radiates in the far infrared (IR) (for a review, see 
e.g., Sanders & Mirabel 1996). For distant galaxies, this emis- 
sion is redshifted to the the submillimeter and millimieter wave- 
lengths. Hence, observations at these wavelengths are regarded 
as key for the studies of the formation and early evolution of 
galaxies. 

The first product of any galaxy survey is the number counts, 
i.e. the number of galaxies as a function of limiting flux in a 
certain sky area. The evolution of galaxies is often investi- 
gated through number counts. Number counts alone do not 
provide non-degenerate information on the redshifts of detected 
sources; therefore one constructs a suitable model and extracts 



the history of galaxy evolution through the model. Various 
models of IR galaxies have been proposed (e.g., Beichman 
& Helou 1991; Franceschini et al. 1994; Pearson & Rowan- 
Robinson 1996; Guiderdoni et al. 1998; Malkan & Stecker 
1998; Tan, Silk, & Balland 1999; Takeuchi et al. 1999, 2001; 
and Pearson 2000) to understand and interpret the observational 
results. Recent studies using the Infrared Space Observatory 
(ISO) (e.g., Kawara et al. 1998; Puget et al. 1999; Dole et al. 
2000; Serjeant et al. 2000; Efstathiou et al. 2000; Kawara et al. 
2000; and Okuda 2000) and the Submillimeter Common-User 
Bolometer Array (SCUBA) (e.g., Smail et al. 1997; Barger et 
al. 1998; Hughes et al. 1998; Barger, Cowie, & Sanders 1999; 
Eales et al. 1999, 2000; Holland et al. 1999; and Blain et al. 
1999, 2000) have stimulated this field. Unfortunately, observa- 
tions in the submillimeter wavebands are difficult, mainly be- 
cause of the following reasons. First, ground-based observa- 
tions in these wavelengths are hard because of the atmospheric 
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absorption. Second, state-of-the-art instruments are required 
both for the antennas and for the detectors. New breakthroughs 
in the submillimeter field are just around the corner with the 
Large Millimeter and Submillimeter Array (LMS A) and the At- 
acama Large Millimeter Array (ALMA). 

The LMSA 6 is the ground-based radio interferometric facil- 
ity proposed by Japan. In the current design concept the array 
will consist of 32 12-m antennas, whose total collecting area is 
3619 m 2 . The baselines can be 20 m to 10000 m, which real- 
izes the maximum angular resolution of 0.0 1 arcsec at 345 GHz. 
The receivers will cover observing frequencies from 80 to 900 
GHz. The array will be located at Pampa la Bola/Llano de Chaj- 
nantor, a very high site in Chile, to realize sub-arcsec resolution 
imaging at the very high frequencies. The Atacama Large Mil- 
limeter Array (ALMA) 7 is the new name for the merger of the 
major millimeter array projects into one global project: the Eu- 
ropean Large Southern Array (LSA), the U.S. Millimeter Array 
(MMA), and the LMSA. Hereafter we use the ALMA to de- 
scribe the whole. The 5er sensitivities of the LMSA at wave- 
lengths 350 /im, 450 /im, 650 /im, 850 /mi, 1.3 mm, and 3.0 
mm (mean values in winter season) are 1200, 660, 370, 49, 23, 
and 14 /iJy beam" 1 , respectively (8-hour integration and 8-GHz 
bandwidth) in dual polarization and DSB mode. For the refer- 
ence of the atmosphere condition and system temperature r sys , 
see Matsushita et al. (1999). The ALMA will consist of 96 12- 
m antennas, and in this case the 5a sensitivities at the above 
wavelengths are 390, 220, 120, 16, 7.5, and 4.6 /dy beam" 1 , 
respectively, under the same assumptions. 

The Atacama Submillimeter Telescope Experiment (ASTE) 
is a Japanese project to install and operate a 10-m submillime- 
ter antenna at Pampa la Bola in the Atacama desert in northern 
Chile. A new 10-m antenna has been constructed in Nobeyama 
Radio Observatory in February 2000, and is scheduled to be in- 
stalled in Chile in 2001. The telescope will be equipped with 
submillimeter-wave SIS mixer receivers and a submillimeter- 
wave direct detector camera. This project has two goals: One 
is to examine the 10-m submillimeter antenna as the technical 
prototype of the LMSA at the site of northern Chile, and the 
other is to explore the possibility of the submillimeter obser- 
vations in the southern hemisphere and to obtain astrophysical 
results. Details of the present development of the ASTE is given 
in e.g., Matsuo et al. (2000) and Ukita et al. (2000) 8 . 

Before the ALMA/LMSA begins operation, a Japanese in- 
frared satellite ASTRO-F will be launched 9 . The ASTRO-F 
will carry out a FIR all-sky survey at two narrow bands, N60 
(50-70 fim) and N 170 (150-200 /im) and two wide bands, 
WIDE-S (50- 110 pm) and WIDE-L (110- 200 /mi). The de- 
tection limits are estimated as 39 mJy and 110 mJy for N60 
and N170, and 16 mJy and 90 mJy for WIDE-S and WIDE- 
L, respectively (Takahashi et al. 2000). A huge database of IR 
galaxies can be expected from the survey. Some other strato- 
spheric and space facilities such as SOFIA (Becklin 1998), 
SIRTF (Fanson et al. 1998), and FIRST (Pilbratt 1998) are also 
in progress, but we stress that the ASTRO-F is the one and only 
all-sky surveyor among these facilities 10 . The combination of 
the FIR data from the ASTRO-F and the submillimeter and mil- 



limeter data from the ASTE and the ALMA/LMSA promises to 
a new era of extragalactic astrophysics and observational cos- 
mology. 

In this paper, we study and examine what we can learn 
about galaxy formation and evolution from the unprecedent- 
edly high-quality data which will be obtained by the ASTE and 
the ALMA/LMSA through a simple empirical galaxy number 
count model proposed by Takeuchi et al. (2001). The rest of 
the paper is organized as follows. In Section 2 we present the 
IR galaxy number count model on which our subsequent dis- 
cussions are based. We also formulate the source confusion 
limit and optimal survey area to suppress both the fluctuation 
caused by the source clustering, and the Poisson shot noise 
by the sparseness of the sources. In Section 3 we show the 
galaxy number counts at submillimeter and millimeter wave- 
lengths and extensively examine how galaxy evolution and for- 
mation epoch affect the number counts. We also discuss the 
optimal strategy of a survey at these wavelength to obtain the 
maximal information on galaxy evolution and formation. Sec- 
tion 4 is devoted to our summary and conclusions. 

2. MODEL DESCRIPTION 
2.1. Number Count Model 

The galaxy number count model is represented by galaxy 
spectral energy distributions (SEDs), luminosity function (LF), 
cosmology, and galaxy evolution. We briefly review the num- 
ber count model by Takeuchi et al. (2001 ; hereafter T01). Based 
on this model, we calculate galaxy number counts with various 
galaxy evolution histories in the subsequent part of this paper. 

The FIR SED is constructed based on the IRAS color- 
luminosity relation (Smith et al. 1987; Soifer & Neugebauer 
1991), 



log |^ = (0.10 ± 0.02) logL 60 - (1.3 ± 0.2) , 

•J 100 



(1) 



where S\ is the detected flux density at wavelength A [/im], 
and L(,q [L Q ] is the intrinsic luminosity evaluated at 60-/im, 
Lgo = vL v . This relation is converted to the dust temperature 
rd ust -L6o relation, then the modified blackbody continuum with 
the corresponding r<j ust are calculated. We then added the mid- 
IR and radio continuum spectrum based on the empirical cor- 
relations between fluxes in these wavelengths and FIR contin- 
uum. The final SEDs that we use in our number count and 
CIRB models are presented in Figure 1 . For the full description 
of the SED construction, see T01. 

The 60-/mi LF based on IRAS data by Soifer et al. (1987) 
is adopted as the local IR LF of galaxies and pure luminosity 
evolution is assumed in this study. Thus the 60-/mi luminosity 
of a galaxy at redshift z is described by 

L 60 (z)=L 60 (0)/(z). (2) 
We also assumed that the luminosity evolution is 'universal', 
i.e. its functional form is fixed. In this paper we call f(z) 'the 
evolutionary factor' as in T01. 

Using the above formulae, we calculate the galaxy number 
counts. As we see later, we set the maximum redshift we con- 
sider the counts of galaxies into the calculation, z max = 5 and we 
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FIG. 1 . — The assumed galaxy spectral energy distribution (SED) in the near infrared to the radio wavelengths. The prominent emission bands are PAH features. 
The SEDs with the FIR luminosity of 10 8 , 10', 10 10 , 10", 10' 2 , 10' 3 , and 10 14 Lq are shown from the bottom in this order. These construction of these SEDs are 
described by Takeuchi et al. (2001). The vertical dotted lines depict the wavebands of IRAS and the two representative wavelengths of the atmosperic window in 
the submillimeter. 



vary it from z max = 2 to z max = 7 to investigate the effect of z max 
on the number counts. T01 statistically estimated the permitted 
range of the evolutionary factor from the observed IR number 
counts and the CIRB spectrum. The evolutionary factor is di- 
rectly connected to the evolution of the cosmic IR luminosity 
density, which is closely related to the cosmic star formation 
history (e.g., Kennicutt 1998). The conversion from the IR lu- 
minosity to the star formation rate and the interpretation of the 
evolutionary factor to the cosmic star formation history are dis- 
cussed in T01 . In this paper we examine how the number counts 
at submillimeter and millimeter wavelengths vary with different 
evolutions, based on the two representative evolutions from T01 
(we call them Evolution 1, and 2 in this paper). 

The model of T01 has been constructed based on the prop- 
erties of IR galaxies, and it slightly underestimates the submil- 
limeter source counts. Therefore in this study, we use an ad- 
ditional model based on the SCUBA source counts. The addi- 
tional evolution model overestimates the CIRB spectrum, but 
the possibility of reconciling the CIRB and submillimeter num- 
ber counts is discussed in Takeuchi et al. (2000); therefore, we 
do not go into the details of this problem further. We call this 
model 'Evolution 3'. These three evolutions are summarized in 
Figure 2. 

2.2. Formulation of Source Confusion 

We summarize the framework of the source confusion statis- 
tics and derive the source confusion limit. The basic formu- 
lation of the source confusion was first discussed by Scheuer 
(1957). The formula for the confusion limit for power-law 
number counts was formulated in the classical work of Condon 
(1974), and extended for the case of general number counts by 
Franceschini et al. (1989). 

In this subsection, we briefly re-formulate these issues and 
examine the formulae carefully for precise numerical calcula- 



tions. The original formulation of Scheuer (1957) was rather 
complex, so we use a simper but modern theory of random sum- 
mation (e.g., Gnedenko & Korolev 1996) to clarify the mathe- 
matical structure of the problem. 

2.2.1. Statistics of Signal from Faint Sources in the Beam 

First we define a differential number count per sr with flux 
Se[S,S + dS],n(S)dS, 

dN 

n(S)dS= — dS . 



dS 



(3) 



where N = N(> S) is the cumulative number count. As a radio 
telescope scan across the sky, the antenna temperature fluctu- 
ates due to many confusing sources passing through the beam. 
Let h(9, (p) be the beam pattern (normalized to unity at the beam 
center), and x = Sh(9, (f>), the response of the telescope to a 
source of flux density S located at an angular position (9, (f>) 
from the beam axis. The mean number of source responses of 
intensity x in a beam is 

R(x)dx = f n (-r^-r) t^— dSl . (4) 

The amplitude of the total signal at any point, D, is a randomly- 
stopped sum of the responses x due to all the sources in the 
response pattern: 

D =x\ +X2~\ Yxk , 

where K is the total number of contributing sources. We should 
note that K randomly varies and takes an integer value. The 
mean number of x's, k, is 

/oo poo 
R(x)dx = / R(x)dx . 
oc JO 

The probability that the total number of x's has a value k is 
given by a Poisson distribution 



K 

Pk = — e 



(5) 
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FIG. 2. — The three representative evolutionary history of galaxies used in this study. The upper and middle models are estimated from the infrared galaxy number 
counts (Takeuchi et al. 2001). The bottom model is the new additional one constructed to reproduce the observed submillimeter galaxy number counts. In the top 
panel, dashed and dotted lines represent the modified evolutionary factors to examine the effect on the evolution at z = 1 ~ 2 to the number count at submillimeter. 



We have the (cumulative) distribution function (d.f.) of x, F(x), 

(6) 



(7) 



F(x) = P(x' < x) = I --—dx' . 

J-oo K 

Then the d.f. of a signal D with k summands is 

r D 

F *(*+D(£)) = / F* k (D-x')dF(x'), 
Jo 

F*\D) = F(D). (8) 
Combining eqs. (5), (6), (7), and (8), we obtain the d.f. of D, 
G(D): 

OO 

G(D) = Y,PkF* k (D). 

The characteristic function (c.f.) of G(D), &c(t), is expressed 
bythec.f. of F(D), $ F (t), as 



/OO ^ 
e i,D dG(D) = Y,Pk^F(ty 



(9) 



k=l 

Recall that pt = (n k e~ K )/k\, then we have 

00 ^k. 



k=l 



k=l 



= e 



k\ 



The expectation value and variance of a randomly-stopped sum 

D are known to be 

E[D] = E[k]E[ Xl ], 

V[Z)] = E[fe]V[x 1 ] + V[^](E[x 1 ]) 2 

(10) 



(e.g., Stuart & Ord 1994) where E[ • ] and V[ • ] represent the 
expectation and variance, respectively. Thus we have 

E[D] = E[k]E[x] = nE[x] 
k I x dx 

K 



-f 

J —o 



xR(x)dx , 



V[Z)] = E[^]V[x] + V[fe](E[x]) 2 
= «(E[(x-E[x]) 2 ] + (E[x]) 2 ) 



/OO 
x 2 R(. 
OO 



x)dx . 



(11) 



(12) 



2.2.2. Confusion limit 

We next formulate the relation between the beam size 9q 
and the source confusion limit. Hereafter we assume that 
h(8, (j>) = h{9). A general case without this assumption is con- 
sidered by Condon (1974), but the result is not substantially 
affected by the above simplification. First we discuss the case 
that the number count is described by a power-law: 



n(S) = aS~'< . 



(13) 



Then we have 
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The mean number of x, R(x) is 



R(x)= I ah(9yx-~>— = an efl x-~> 



where Sl e ff is 



eff = / h(6)~<- x dtt. 



(15) 



(16) 



Here we consider the Gaussian beam pattern: 

h(6) = e- (4ln2 <%Y . (17) 

We note that more complex beam patterns can be taken into 
account. Other beam patterns have been commented upon by 
e.g., Scheuer (1957) and Franceschini et al. (1989). But in a 
real situation, the result does not depend on the details of the 
beam pattern, and a Gaussian beam can be used as a represen- 
tative case. In Equation (17), 9q is the FWHM of the beam and 
relates to the standard deviation of the Gaussian, a, through 
assume that the beam area on the sky is 
small enough that the solid angle can be integrated on a plane. 
Then Condon (1974) derived the effective beam size as 



^eff = 



2tt 



-(41n2)( 7 -l) 



9d9d(p 



I 



o 

oo 



-(41n2)(7-l)i/> 



dip 



(41n2)( 7 -l)' (18) 

(1 < 7 < 3). In the third step we set tp = (9/9 ) 2 . We obtain the 
confusion limit flux to a cutoff deflection x c : 



cr(x c ) 



f 

JO 



x R(x)dx 



aQ e ffX 2 1 dx = afl e ff / x 2 1 dx 







and thus 



O^eff 

3^ 



a(x c ) = 



,3-7 



eff 



3-7 



(19) 



(20) 



As usually defined, if we set x c = qa, we have 

r i 



3-7 



3-T 3--y 

q 2 a 2 



Thus we obtain the important expression for the case of a 
power-law number count, 



3-7 



For the Gaussian beam pattern, eqs. (21) and (18) lead to 



3-7 



(21) 



(22) 



^3-7/ V( 41n2 )(7-1)/ 
This formula is used very frequently, probably because it is ex- 
pressed in a simple analytic function. However, we should note 
that this derivation includes a mathematical flaw because 9 is 
integrated from to 00, even though 9 must be small. Fortu- 
nately, the above integration converges, but when we general- 
ize the form of the number counts, it diverges, as we discuss 



below. Therefore we should set a certain cutoff in the integra- 
tion in the real calculation. Thus, the confusion limit given by 
Equation (22) is too pessimistic. 

We now turn to formulate the confusion limit for the general 
number counts (Franceschini et al. 1989). Equation (4) leads 
to 

n( X \ dn ~ 



x R(x)dx- 



I * 2 *J. "( 



-(41n2)V> 



\h(9) J h{9) 

dip 



dx 



-(4ln2)i/> 



dx 



4ln2 



n{rjx)drj 



dx : 



4ln2 



Kxc) , (23) 



where ip c and r) c are the upper cutoff, which corresponds to the 
beam area and is set to obtain a numerically reasonable result. 
The result strongly depends on the value of the cutoff. Just as 
in the above discussion, we set x c = qa, so 



irHqa) 



{) 



(24) 



4ln2 

Thus, we obtain a general relation between the beam size and 
the confusion limit as 



6*o = 



/ 4ln2 
Trl(qcr) 



(25) 



2.3. Required Survey Area 



Next we formulate the calculation for the necessary sky area 
of the survey for the studies of galaxy evolution. In order to 
estimate the effects of galaxy evolution from survey data, a sig- 
nificant sky area should be scanned to suppress the variation in 
the galaxy surface number density. Since galaxies are clustered 
on the sky, the nominal error bar estimated from the assump- 
tion of Poisson statistics is an underestimate, and we must take 
the galaxy angular correlation function into account. Barcons 
(1992) considered the effect of clustering in his fluctuation anal- 
ysis. Despite its importance, this problem is often overlooked 
in the literature on deep surveys. At the faintest fluxes, this ad- 
ditional effect is gradually diluted by the projection along the 
line of sight. We formulate the effect in the following. 

Consider a survey area 51 and divide it into small cells { AO,} 
so that the number of galaxies in the cell {AO,}, Mi = or 1. 
We set the mean galaxy surface number density Af. Then we 
have a mean number in a solid-angle cell Afi„ Af if 

(M) = TVAfij . (26) 
By definition, we have (A//) = (Af) = (A^ 3 ) = • • •. We observe 
(MiJVj) = A^( 1 + w(9ij)) AQ, AO; (27) 
((M- (M»(A0- (A/}») = (MA/}) - (AO (A}) 

= Af 2 w(9 ij )An i ACl j , (28) 
where w(9) is the angular two-point correlation function of 
galaxies. Next we consider the number of galaxies N in the 
survey area We have 



(N)= (M>= / A"d^ = Al7 



(29) 



i j 

= [ Afdn+ [ [ AA 2 (l+w(6>i2))dQidn 2 
Jn J Jn 



■■Afn+Af 2 n 2 +M 2 



w(9 12 )dn 1 dn 2 . (30) 
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and thus 



((N-(N)) 2 ) = (N 2 )-(N) 2 =Arn + Af 2 [ [ w(6 12 )dSlidQ&l) 

J Jn 

If we assume that we take sufficiently wide area whose dimen- 
sion is much larger than the coherent scale of angular cluster- 
ing, we can approximate the above expression as 

{(N-{N)f) ~NVL+N 2 n ( w(6)dQ. 

Jn 

= Afn(l+AfJ w(6)dnj . (32) 

Here we evaluate the 'signal-to-noise ratio' of the number 
counts S/N: 



S/N = 



y/{(N- (n)) 2 ) ^jvn(i+Nf n w(6)dn) 



J n w(6)dQ ' 



(33) 



The approximation used in the third step can be adopted when 
A/T2 is sufficiently large. We see that the S/N depends on the 
angular correlation strength and the solid angle of the survey 
and is almost independent of the surface density of the sources 
W. This fact shows that, in order to determine galaxy evolu- 
tion from number counts, large-area surveys are required. The 
angular correlation function of the IRAS galaxies is w(9) = 
(6/6 y 0M (6*o = 0? 1 1 : Lahav, Nemiroff, & Piran 1990). We as- 
sume that the submillimeter sources and the faint radio sources 
have the same clustering properties as IRAS galaxies. But even 
if we observe the same species of galaxies as the IRAS galax- 
ies, an IR galaxy appears to be bright at FIR but faint at submil- 
limeter; therefore, we should assume their SEDs and convert 
the correlation function. For the relatively bright radio sources, 
a steeper clustering power-law exponent w(8) oc #-( 11_L4 ) has 
been reported by Cress et al. (1996). The angular correlation 
of galaxies out to a fainter flux limit is obtained by using the 
scaling relation of w(8) with a detection limit through the rela- 
tivistic Limber's equation (e.g., Peebles 1980). 

Using this formula, in Section 3.3 we estimate the minimum 
survey area required for a sufficient S/N. 

3. RESULTS AND DISCUSSIONS 

3.1. The Galaxy Number Count Predictions Expected in the 
ALMA/LMSA 

First we show the galaxy number counts [sr _1 ] from 90 fim 
to 3 mm expected by the ALMA/LMSA in Figure 3. The dot- 
ted curves describes the no-evolution predictions at each wave- 
length. The other curves depict the expected number counts 
with the three evolutionary histories discussed in Section 2.1, 
Evolution 1, 2, and 3. All of these evolutions give a good fit 
to the FIR data and we cannot distinguish between each other 
with the present quality of the observed data (T01). Some ob- 
served source counts are also shown in Figure 3. In the IR, 
we plot the galaxy counts from Stickel et al. (1998), Kawara 
et al. (1998), Puget et al. (1999), and Dole et al. (2000), Juvela, 
Mattila, Lemke (2000), Matsuhara et al. (2000), and Efstathiou 
et al. (2000). We also show the IRAS 100 /im galaxy counts by 
a hatched thin area at brightest flux. In the submillimeter, the 
data are taken from Smail et al. (1997), Barger et al. (1998), 
Barger, Cowie, & Sanders (1999), Hughes et al. (1998), Eales 



et al. (1999), Holland et al. (1999), and Blain et al. (1999, 
2000). At 1 .3 mm, we show the preliminary result of MAMBO 
(Max-Planck Millimeter Bolometer array) reported by Bertoldi 
(2000). 

3.2. The Source Confusion Limits 

We calculated the source confusion limits at the submillime- 
ter and radio wavelengths based on Evolution 2. The results are 
shown in Figure 4. The angular resolution of the ALMA will be 
better than 0. 01, i.e., there will never be a confusion problem 
in the ALMA observation, and we do not show it in Figure 4. 
When the slope of the number counts is flat, the source con- 
fusion limit will be much improved if better angular resolution 
can be achieved. On the other hand, in the case that the slope of 
the number counts is very steep, the source confusion limit will 
not be drastically improved even with a much better angular 
resolution. Since the slope of the submillimeter counts is steep 
in the brighter flux regime and gets flatter towards the fainter 
fluxes, the confusion limit becomes lower quickly if the beam 
size (FWHM) becomes small. 

In the submillimeter regime, we see that SCUBA 
(beam size = 14.5 arcsec) deep surveys are nearly confusion- 
limited. Recently Hogg (2000) pointed out that if the number 
count slope is steep, i.e., dlogN / dlogS <^ -1.5, the informa- 
tion of faint sources can be completely destroyed by the confu- 
sion effect. Eales et al. (2000) also performed a Monte Carlo 
simulation and concluded that their estimated fluxes of the sub- 
millimeter SCUBA sources are actually boosted upwards sig- 
nificantly. The median boost factor is 1 .44, with a large scatter; 
in the worst case the bias reaches an order of magnitude. They 
also showed from their simulations that 19 % of their SCUBA 
sources have positional errors greater than 6 arcsec. Therefore, 
we should be quite cautious of the fluxes or positions of the 
faintest submillimeter sources. Considering this point, we can 
obtain the flux safely for sources brighter than ^10 mJy. In the 
next subsection, we see that a reliable survey at brighter fluxes 
is crucially important for the study of galaxy evolution. The 
large-area submillimeter survey of sources brighter than 10 mJy 
by the ASTE is suitable for this purpose. 

On the other hand, the high angular resolution of the 
ALMA/LMSA results in a very low confusion limit of < 1 /iJy. 
Thus, we can estimate the flux and position of such faint sources 
very precisely. Since the slope of the number counts of millime- 
ter sources becomes flatter at such faint fluxes, improving the 
angular resolution provides a drastically better confusion limit. 

3.3. Submillimeter Bright Source Counts and Galaxy 
Evolution at z = 1 ~ 2 

In order to derive information on evolution at z = 1 ~ 2 from 
submillimeter source counts, we modify the evolutionary fac- 
tor and examine how the number counts vary with the modifi- 
cation. We show the evolutionary factors used here in the top 
panel of Figure 2. The thick solid step function is Evolution 1 
itself. We first modify Evolution 1 so that the evolutionary fac- 
tor at z = 1 - 2 is 10 1 5 . Secondly we set the evolutionary factor 
at z = 1 -2 to be 10 2 . The resulting number counts are shown 
in Figure 5. The effect of the different evolutions at z = 1 - 2 is 
clearly seen in the bright end of the number counts in Figure 5. 
However an important problem arises here. For measuring the 
number counts in such a bright flux regime, a very large-area 
sky survey is required. Since the number density of the sources 
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FIG. 3. — The galaxy number count predictions from the infrared to the radio wavelengths. The dotted curves represent the number counts of galaxies without 
evolution. The dot-dashed, solid, and long-dashed curves describe the number counts with Evolution 1, 2, and 3, respectively. In the infrared (90 and 170 fim), we 
plot the galaxy counts from Stickel et al. (1998), Kawara et al. (1998), Puget et al. (1999), Dole et al. (2000), Juvela, Mattila, & Lemke (2000), Matsuhara et al. 
(2000), and Efstathiou et al. (2000). We also show the IRAS 100 fim galaxy counts by a hatched thin area at brightest fluxes. In the submillimeter (350, 450, 650, 
and 850 /im), the data are taken from Smail et al. (1997), Barger et al. (1998), Barger, Cowie, & Sanders (1999), Hughes et al. (1998), Eales et al. (1999), Holland 
et al. (1999), and Blain et al. (1999, 2000). At 1.3 mm, the preliminary result of MAMBO (Max-Planck Millimeter Bolometer array) reported by Bertoldi (2000) is 
plotted. The vertical dot-dot-dot-dashed lines represent the 5<7-detection limits of the LMSA, and vertical dashed lines shows the 5<7-detection limits of the ALMA. 



is so small and the statistical fluctuations are large, secure de- 
termination of the counts is difficult. 

The necessary sky area is estimated by using the formulation 
presented in Section 2.3. The characteristic depth of a survey is 
defined by the survey flux limit, S um . As we mentioned in Sec- 
tion 2.3, we here assume that the submillimeter sources brighter 
than 100 mJy are clustered as strongly as the IRAS sources 
brighter than 0.7 Jy, which have the angular correlation func- 
tion w Q (9) = (9/9qT 0M . 

We define the relative characteristic depth of a survey by 
d* = (0.1 [Jy]/5ii m [Jy]) 1 / 2 . Then, the angular clustering w(9) 
is scaled as 



w(6) = d; 1 wo(d*6) 



(34) 



where wq(9) is the angular correlation function of the survey 
whose flux limit is 0. 1 Jy. Under the assumptions mentioned 
above, the clustering of the submillimeter sources is expressed 
as 

w(0) = d: lM w Q (9). (35) 



The signal-to-noise ratio defined in Section 2.3 becomes 



S/N ' \<i. ( -".L^<> 

and we obtain 



Slim = 0.1 



05 '-(.V:„, 0.1)^T° 



f n w (6)dn 



1/0.83 



(36) 



^S/N) 2 J nW0 (9)dn / 

In the left panel of Figure 6 we show the area-limiting flux 
relation derived directly from eqs. (33) and (35). In fact, the 
equation is well approximated by 

ft(S lim ) cx 5f im . (37) 

The above approximation is valid as long as the fluctuation by 
clustering dominates the variance. Since recent observations 
of the submillimeter number counts suggest strong galaxy evo- 
lution (Section 1), the number counts increase rapidly toward 
fainter fluxes. Consequently, the Poisson fluctuation is small 
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FIG. 3. — Continued. 



as compared with the clustering variance and the condition for 
the approximation is satisfied. Thus, the relation presented in 
Figure 6 is well represented by Equation (36). 

We can also estimate the relative observation time to per- 
form such a survey by multiplying S, 2 im by the obtained sky area 
£!(Sii m ). It is straightforward to see that the relative observation 
time is roughly proportional to V-Siim- Therefore, the fainter the 
limiting flux, the shorter the required observation time. This 
is presented in the right panel of Figure 6. But actually, as 
we discussed above, we cannot make the flux limit too faint 
(<; 10 mJy) because of severe source confusion. 

If we adopt S/N = 5, the required survey area is 10 deg 2 
at 10 mJy, which is large but not unreasonable for a present- 
day submillimeter survey. If S/N = 3, then the survey area is 
20 deg 2 at 50 mJy. Thus, the large-area survey by the ASTE 
is suitable for this purpose. For the minimal survey of 1 deg 2 , 
the virtual observation time is estimated to be approximately 
200 hours by a single-pixel detector. If the array-type detector 
is available, this survey will be performed much more quickly 
and efficiently. 

We note that these predictions depend on the assumed clus- 
tering strength. When the survey is finished, we will have a 
large database of submillimeter sources with firm counterpart 



associations. Using the resulting catalog, we will be able to 
estimate the angular clustering properties of the submillimeter 
sources. If the variance of the detected source counts are differ- 
ent than what we assumed, then it provides important quantita- 
tive information on the clustering evolution of dusty galaxies. 
We stress that this will be another important purpose of this 
survey. 

3.4. Faintest Source Counts and Galaxy Formation 

We next adopt three other cases for the evolutionary factor 
f(z) to examine the effect of the galaxy formation epoch on the 
number counts at submillimeter and radio wavelengths. First 
we use the evolutionary histories based on Evolution 2, and 
vary the the cutoff redshifts as Zf orm = 2, 3, 5, and 7. We note 
that this cutoff redshift represents the epoch of the formation or 
appearance of the dusty galaxies, and does not necessarily im- 
ply the formation of all species of galaxies. In addition, we also 
study the case of submillimeter-based Evolution 3, introduced 
in Section 2. 1 . For the modified evolutionary histories based on 
the Evolution 3, the redshift cutoff are also set to be Zf orm = 2, 
3, 5, and 7. 

The resulting number counts with various redshift cutoffs are 
shown in Figure 7. The top two panels present the number 
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FIG. 4. — The 5cr-source confusion limits as a function of the beam size, calculated from the Evolution 2 in Figure 2. The solid, dotted, dashed, and dot-dashed 
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FIG. 5. — The resulting number counts produced by the evolutionary factor in the top panel of Figure 2. The dotted line shows the no-evolution result. The solid 
curve represents the galaxy number counts with Evolution 1 . The dashed and dot-dashed lines represent the number counts with the modified evolutions so that the 
evolutionary factor at z = 1 -2 is 10 1 ' 5 and 10 2 , respectively. The effect of the evolution at z = 1 -2 is clearly seen in the bright end of the number count. 



counts at 850 ^iin and 1 .3 mm with Evolution 2. The bottom 
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FIG. 6. — Left panel: The required sky area to suppress the galaxy clustering fluctuation as a function of limiting flux of a survey. Right panel: Relative observing 
time to perform the corresponding survey. The unit of the time is arbitrary. 



tow panels show the number counts with Evolution 3. The ef- 
fect of the cutoff is clearly seen in the faintest end of the number 
counts in these two panels. In the bottom panel of Figure 7, the 
effect of redshift cutoff is more prominent than in the top panel. 

In contrast to the case discussed in Section 3.3, a very deep 
detection limit is required to investigate the redshift cutoff, i.e., 
the galaxy formation epoch. The extremely deep survey by 
ALMA/LMS A is a unique opportunity for this purpose. At such 
a faint flux level, the variance of the number count is diluted by 
the large source surface density and the source projection; we 
can, therefore, distinguish between the predictions of different 
models presented in Figure 7. 

When we discuss such a deep flux limit, the dynamic range of 
a detector should be taken into account. If a very bright source 
exists in the field of view, the faintest sources near the detection 
limit cannot be detected in the field. We evaluated the proba- 
bility that the sources 10 3 times brighter than the 5a-detection 
limit of the LMSA and ALMA exist in a field of view. The 
field of view of the LMSA observation is 20-arcsec diameter 
circular area (7.4 x 10~ 9 sr). The number density of the sources 
10 3 times brighter than the 5<r limit at each waveband, N, is 
directly obtained from Figure 3. For example, N of the ALMA 
at 850 ^m is 4 x 10 4 sr" 1 . Then we obtain the expectation value 
of the number of very bright sources compared with the detec- 
tion limit found in the field of view, /j, = N x (field of view). By 
using this value, we can treat the probability of finding k bright 
sources in the field as a Poisson process 

W=-fcf- ■ (38) 

Then the probability that more than one source exist in the field 
of view is 

1-^(0)= 1-e-". (39) 

Even in the worst case of the ALMA (when we use Evolution 3) 
the probability is <~ 3 x 10~ 4 . Thus, the probability that a very 
bright radio source lies in the field of view is so small that we do 
not have to worry about the dynamic range, and we can safely 
obtain the information of very distant faint sources. 

3.5. Redshift Estimation from the Dust Continuum 



Determination of source redshifts is an important but quite 
challenging issue in cosmological studies. Information on 
source redshifts enables us to construct the LF as a function 
of redshift. We here assume the LF shape to be unchanged 
during evolution, but this is a strong assumption which should 
be tested with observational data. Spectroscopic observations 
will be extremely difficult or almost impossible on very dis- 
tant objects, and alternative methods are necessary. In the op- 
tical wavelength, photometric redshift technique has provided 
a breakthrough for further studies of ultra high-redshift objects 
(e.g., Fernandez-Soto et al. 1999). 

On the other hand, at longer wavelengths such as the far-IR 
(FIR) or submillimeter it is harder to estimate source redshifts 
because of the smooth nature of their SEDs. Takeuchi et al. 
(1999) tried to estimate redshifts roughly from FIR photometry 
and found a practical possibility. Carilli & Yun (1999) proposed 
the radio-to-submillimeter spectral index as a redshift indicator 
for star-forming galaxies. They used the correlation between 
the FIR flux from the thermal dust and the radio flux from the 
synchrotron emission from the interstellar matter and supernova 
remnants, and utilized the spectral index between the two kinds 
of emission as a function of redshift. 

Though the radio-to-submillimeter index method is now of- 
ten used to estimate the redshift of submillimeter sources, the 
origin of the two components of radiation is substantially dif- 
ferent. If we can estimate the redshift of the sources from the 
ratio of two or more flux densities which have the same physical 
origin, it will be undoubtedly the most desirable method. 

We here develop such a method which is an extension of 
the attempt of Takeuchi et al. (1999). We use FIR and sub- 
millimeter/millimeter flux densities for this purpose. For the 
FIR database, we suppose the wideband photometric catalog 
of the ASTRO-F FIR all-sky survey. At FIR to submillimeter 
wavelengths, continuum emission is dominated by blackbody 
radiation from the big dust grains in equilibrium with the am- 
bient ultraviolet radiation field. Thus, the flux densities in this 
wavelength range are radiated from the same emission mecha- 
nism. 

First we show the redshifted SED of an infrared galaxy 
with luminosities 10 8 , 10 10 , 10 10 , and 1O 12 L at redshift z = 
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FIG. 7. — Top panels: The galaxy number counts at 850 (im and 1.3 mm with redshift cutoffs Zf orm = 2,3,5, and 7 based on Evolution 2 in Figure 2. Bottom panels: 
The same as the Top panels except that the number counts are based on Evolution 3 in Figure 2. In both panels, the dot-dashed, dashed, solid, and long-dashed 
curves represent the resulting number counts with redshift cutoffs zf orm = 2,3,5, and 7, respectively. A very low detection limit is required to investigate the redshift 
cutoff, i.e., the galaxy formation epoch. 



0.1,0.5,1.0, and 5.0 in Figure 8. The 5cr-detection limits of 
ASTRO-F and LMSA are also shown. The problem is that the 
peak of the blackbody radiation shifts not only with redshift but 
also with dust temperature r aust . Therefore, this method suffers 
from this degeneracy. Is it useless as a redshift estimator? In 
order to clarify this point, we present the detected flux densities 
and flux density ratios of star-forming galaxies as a function of 
redshift in Figure 9. The top two panels are the flux-redshift re- 
lations at 170 fim and 850 /im. The horizontal dotted line in the 
top-left and top-right panels depicts the detection limit at each 
wavelength. The middle and bottom panels present the color- 
redshift relations. Submillimeter colors are a strong function of 
redshift, and we expect that they work well as rough redshift 
indicators. 

But the ambiguity is of order unity and is too large for, e.g., 
estimating the shape evolution of the source LF. Here we can 
use some additional information. We use the empirical re- 
lation between FIR luminosities of galaxies and the flux ra- 
tio S(,Qfim/Sio()tim (Equation (1)) when we construct the SED 
model. Therefore when we obtain the flux ratio of a source, we 
can estimate the luminosity L v (at 60/mi). Then we compare the 
estimated flux density 



m ■ 



(l+z)L„(at60 /im) 
47rc/ L (z) 2 



(40) 



with the detected flux density at the corresponding wavelength. 
If the two flux densities are significantly different, we correct 
the assumed r<j ust and repeat this trial. Thus we can obtain 
the redshift estimation by an iterative procedure. Galaxies do 
not degenerate in the color-color-flux three-dimensional space 
(Figure 10), and we safely use this method. We call this the 
dust-z method. 

In order to examine the accuracy of the dust-z method, we 
performed a series of Monte Carlo simulations. We assumed 
the following observational wavebands: 60, 90, 170, 450, 850, 
1300, and 3000 /im. The errors in flux measurements are set to 
be 5 % and 30 %. We input the assumed SEDs and added ran- 
dom errors to them, and calculate their 'observed' fluxes. Then 
we performed the above algorighm, and estimated the redshifts 
and luminosities of simulated galaxies. The result of the sim- 
ulation is presented in Figure 1 1 . We show the estimates for 
galaxies with input luminosity 10 10 , 10 12 , and 10 14 L Q . The in- 
put redshifts of galaxies are fixed to 0.5, 1.0, 2.0, and 5.0. We 
examined two representative cases where the errors are ~ 5 % 
and <~ 30 %. In these simulations, we calculated 100 realiza- 
tions for each redshift and luminosity. We conclude that if the 
error in each band is ~ 5 % (i.e. S/N <~ 20), then the redshift 
can be successfully estimated by the dust-z method. Even in 
case the error is ~ 30 %, the uncertainty in the estimation is 
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comparable or even better than some other presently used indi- 
rect methods for redshift estimation (cf. Rengarajan & Takeuchi 
2001). 

We also evaluate the uncertainty in the adopted SEDs. Tak- 
ing the error in Equation (1) into account, we assess the error of 
~ 40 % in luminosity estimation. The luminosity of galaxies in 
the FIR is well approximated by Lfir oc hence the error 
in temperature estimation is at most (1 iO.4) 1 / 4 " 1 " 7 ~ 1 ± 10 %. 
The peak frequency, ^ pea k, obeys the Wien displacement law: 
^peak = 5-88 x 10 10 T. Thus, the error associated with the uncer- 
tainty in SED templates is evaluated as ~ 10 %. 

This method works effectively when we have many photo- 
metric bands or channels. Therefore, instruments that have a 
large number of wavebands are very useful for this purpose. 
This method provides us not only the redshift information but 
also the dust temperature of the submillimeter sources at the 
same time. This will surely be a strong constraint on the star 
formation history of galaxies. 



3.6. 1 A-GHz Radio Source Counts 

We show the comparison of the observed 8, 5, and 1.4 GHz 
source counts and our model counts in Figure 12 to examine the 
contribution of star forming galaxies to the faint radio galaxies. 
The observational data are taken from Windhorst et al. (1995), 
Richards et al. (1998) (8 GHz), Fomalont et al. (1991), Hammer 
et al. (1995) (5 GHz), White et al. (1997), Ciliegi et al. (1999) 
and Gruppioni et al. (1999) (1.4 GHz). At 8 GHz the observed 
source counts are slightly higher than our model prediction, but 
the discrepancy is not so significant because the survey areas 
are quite small in these studies. The 5 GHz data are well re- 
produced by the present model. Recent studies revealed that 
the faint radio sources with flux densities S v ~ 1 - 10 /iJy are 
dominated by actively star-forming galaxies (e.g., Haarsma et 
al. 2000). At 1-0.01 Jy we see an excess of 1.4-GHz sources 
compared with our model calculations. The excess component 
mainly consists of radio galaxies that are dominated by ellipti- 
cals and are not directly related to the star formation in galaxies. 
But the observed counts and model predictions show a conver- 
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FIG. 9. — The flux-redshift and color-redshift diagrams of the model infrared galaxies with luminosities same as Figure 8. Top two panels are the flux-redshift 
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panels present the color-redshift relations. The vertical dotted lines represent the redshift z = 1 and 2. 



gence with each other toward the fainter flux regime. This sup- 
ports the claim of Haarsma et al. (2000); thus, we expect that 
the millijansky radio sources are really star-forming galaxies. 

In such a long wavelength as 1 .4 GHz (20 cm), we must be 
careful of the diffraction limit. We show the 5cr-confusion limit 
of the 1.4 GHz observations in Figure 13. The model basically 
treats the contribution of the star-forming galaxies, but as we 
mentioned above, bright nonthermal sources contribute to the 
source confusion at this wavelength. We take their contribution 
to the confusion estimation. The dashed line in Figure 13 is the 
confusion when only the star-forming galaxies are taken into 
account, and the solid thick line represents the confusion limit 
including strong nonthermal radio sources. 

Then we investigate what we know from the faint radio 
source counts. We present the radio number counts based on 
the evolutionary histories presented in Sections 3.3 and 3.4 in 
Figures 14 and 15. Figure 14 clearly shows that the faint radio 
counts at 1 .4 GHz depend on the evolutionary status of galax- 
ies at z = 1-2. On the contrary, Figure 15 demonstrates that 
the 1.4-GHz counts are almost insensitive to the evolutionary 



status of galaxies at z > 2. The redshift information by direct 
measurement is required for very high-z objects to evaluate the 
star formation activity, as studied by Haarsma et al. (2000). 

This result shows that deep radio surveys are another probe 
of galaxy evolution in the redshift range z = 1-2. It should 
be tested by comparison with future results from submillime- 
ter large-area surveys, discussed in Section 3.3. 

4. SUMMARY AND CONCLUSIONS 

In this paper we investigated what we can learn about galaxy 
formation and evolution from the data which will be obtained 
by the forthcoming new large submillimeter/radio facilities, 
mainly by the ASTE and the ALMA/LMSA. 

We first calculated the source counts from 90 fim to 3 mm 
by using the infrared galaxy number count model of Takeuchi 
et al. (2001). Based on their number counts, we then derived 
the source confusion noise and estimated the confusion limits 
at various wavebands as a function of the characteristic beam 
size. 

We found that, at the submillimeter wavelengths, source con- 
fusion for the 10 - 15-m class facilities becomes severe at the 
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FIG. 10. — The color-flux diagrams of the redshifted model infrared galaxies at various luminosity and redshift. 
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10 14 Lq, respectively. The vertical dotted lines depict redshifts of 0.5, 1.0, 2.0, and 5.0. Left panel shows that the added random errors are ~ 5 %, and right panel 
presents the case that the errors are ~ 30 %. 



0.1 to 1 mJy level, and astrometry and flux measurement are 
difficult. Then we showed that a very large-area survey of the 
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FIG. 12. — The comparison of the observed 8, 5, and 1.4 GHz source counts and our model counts. The same as Figure 3, the dotted curves indicate the number 
counts of galaxies without evolution, and the dot-dashed, solid, and long-dashed curves represent the number counts with Evolution 1, 2, and 3, respectively. The 
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submillimeter sources brighter than 10-50 mJy can provide a 
unique constraint on the infrared galaxy evolution at z = 1-2, 
and such a survey is suitable for a facility such as the ASTE. 
A large area is required to suppress the statistical fluctuation 
caused by galaxy clustering on the sky. Such a survey will also 
enable us to study the clustering properties of the bright sub- 
millimeter sources, which is still highly unknown. 

We also found that the 5cr-confusion limit of LMSA reaches 
to 1 /iJy, which enables us to study the contribution of sources 
at extremely large redshifts. The source counts at such a faint 
flux level give important information on the galaxy formation 
epoch. At such faint fluxes the statistical uncertainty stated 
above is small enough that we can safely estimate the galaxy 
number counts. 

We then discussed the possibility of using multiband pho- 
tometric measurements in the infrared (by ASTRO-F in this 
work) to the millimeter as a rough redshift estimator. We sug- 
gested that the source redshift can be obtained from the color 
and flux of the thermal dust emission. More precise redshift 
values can be estimated by 'the dust-z method', which addi- 
tionally uses the relation between the Seo^m/Sm^m c °l° r an d 
the luminosity at 60 ^m. Thus multiband instrument is very 



welcome for this purpose. We examined the effectiveness of 
this method by Monte Carlo simulations and found that it suc- 
cessfully works if we have reasonable measurement accuracy. 
In addition, the forthcoming IR facilities such as ASTRO-F and 
SIRTF will provide very good information of IR SEDs of galax- 
ies, which will certainly improve the validity of the method. 

We calculated the comparison of the observed 1.4 GHz 
source counts and our model counts, to examine the contribu- 
tion of star forming galaxies to the faint radio galaxies. The 
faint radio counts at 1 .4 GHz depend on the evolutionary status 
of galaxies at z = 1-2 but are insensitive to the evolutionary sta- 
tus of galaxies at z > 2. In order to explore the radio properties 
of such high-z sources, we need a direct measurement of their 
redshifts. 

First we thank the anonymous referee for useful suggestions 
and comments. We also thank Drs. Hiroshi Shibai, T. N. Ren- 
garajan, Hiroshi Matsuo, Hideo Matsuhara, Nobuharu Ukita, 
and Tomonori Totani for fruitful discussions and comments. 
We are grateful to Drs. Shin Mineshige and Yasushi Suto for 
continuous encouragement. HH, KY, and KN are grateful to 
JSPS fellowship. We made extensive use of the NASA's Astro- 
physics Data System Abstract Service (ADS). 
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FIG. 13. — The 5cr-source confusion limits as a function of the beam size at 1.4 GHz. The same as Figure 4, we calculated the limit from the Evolution 2 in 
Figure 3. The contribution of nonthermal radio galaxies at the bright flux regime is also involved in the calculation. 
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